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THE MASS distribution OF BLACK HOLES IN BL LACERTAE OBJECTS 
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ABSTRACT 

In this Letter, we make an attempt to estimate the masses of black holes and their distributions in BL Lac 
objects. We find there is a bimodal population of black hole masses in BL Lacs: the first has a lower-mass peak 
centred at ^ 7.0 x lOf'MQ composed of high frequency-peaked BL Lacs (HBLs) and the second is a higher-mass 
(> K^Mq) population composed of low frequency-peaked BL Lacs (LBLs). The low mass family of black holes 
in HBLs does not satisfy the Magorrian et al. relation whereas the high mass family in LBLs does this relation and 
is in agreement with Laor's limit. The origination of the bimodal distribution of black hole mass remains open. 

Subject headings: galaxy: active - BL Lac object - black hole: mass distribution 



1. INTRODUCTION 

The black hole masses and their distributions in active and in- 
active galaxies are of great interest in astrophysics (Richstone 
et al. 1998, Kormendy & Gebhardt 2001, Merritt & Ferrarese 
2001). BL Lacertae objects are characterized by their powerful 
featureless continuum radiation spanning from radio to gamma- 
ray bands (Bregman 1990). They play an important role in the 
evolution and unification scheme of active galactic nuclei (Urry 
& Padovani 1995, Antonucci 1993). With extensive studies of 
the BL Lacs, there has been considerable progress in under- 
standing the properties of their multiwaveband continuum and 
radiation mechanisms. It is generally believed that synchrotron 
self-Compton (SSC) emission of high energy electrons in a rel- 
ativistic jet is responsible for the multiwaveband continuum of 
the BL Lacs (Urry & Padovani 1995). The emission is highly 
beamed. The featureless continuum of the BL Lacs only per- 
mits us to use timescale of variability Af to estimate black hole 
mass, Mbh < c^VAt/G, where c is the speed of light, G the 
gravitational constant and V Doppler factor However this only 
provides a very crude upper limit of a black hole mass since it 
is very difficult to determine the shortest timescale, especially 
when there is an uncertainty of relativistic boosting factor. The 
host galaxies of BL Lacs have been discovered by HST, show- 
ing they are normal giant ellipticals (Scarpa et al 2000, Urry et 
al. 2000). Urry et al. (2000) argue that the host galaxy can not 
determine the central engine of BL Lac objects because there 
is no differences between HBLs and LBLs. There is no mea- 
surements of black hole masses in BL Lac objects from stellar 
dynamical method. 

It has been pointed out that the masses of black holes in all 
the radio-loud quasars of low redshift PG sample should exceed 
IO^Mq (Laor 2000). This suggests the black hole mass con- 
trols the radio-loudness R, defined as the ratio of radio to opti- 
cal fluxes. The relativistic jet is generally thought to originate 
from the innermost vicinity of black hole (Blandford 2001). 
The emission from the jet may reflect the black hole parame- 
ters in a more direct way. There is a large sample of BL Lacs 
that the peak frequencies and luminosities are measured. This 
allows us to study the masses of black holes in BL Lacs. 

In this Letter, we present a way to link the black hole mass 
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to the peak frequency and luminosity in BL Lac objects. This 
method allows us to statistically study the mass distribution of 
black holes in a large sample. We find there is a bimodal popu- 
lation of black hole masses in BL Lac objects, corresponding to 
HBLs and LBLs, respectively. The implications of this bimodal 
distribution and Magorrian et al. relation are discussed. 

2. LIMITS ON SUPERMASSIVE BLACK HOLE MASSES 

We intend to connect the mass of central black hole with 
the two available observables, the peak frequency and its 
luminosity Lpk, rather than to fit the continuum. As the first 
step, we express the Doppler factor as a function of i/p^ and 
Lpk. For the simplicity, we assume a spherical geometry for the 
radiation region and a monoenergetic distribution of relativis- 
tic electrons responsible for the continuum emission, namely, 
n(7) = no<5(7-7o), where no is the electron density and 70 the 
Lorentz factor of electrons. Using the standard formula for 
synchrotron emission in Pacholczyk (1970), e^i = noc^B^Fix), 
where C3 = 1 .87 x 10"^^, we have the peak luminosity Lpk in an 
observer's frame 

L,i, = V'^7rRX^e,,^, (1) 

where J-(x) = x Ks/2(z)dz, ^5/3(2) is the modified Bessel 
function with of order 5/3, x = v' lv[ and v[ = 4.2 x 10*B7o. 
The primed parameters are in the co-moving frame of the jet. 
The peak frequency is given by d [x!F(x)] /dx = 0, namely, 

1 /2 

[2 J^(.Xpk)/ /3 (jCpk)] , we have 



z/pk = Vv'^^ = Pxpkz^; = 5.5 X 10''PB7o, (2) 



where jCpk « 1.32. The radiating region must be optically thin 
for synchrotron self-absorption at the peak frequency, namely, 

rsyn(j^') = RoC4B^^^nQiy'~^^^Ks/2(xp^) = tq, (3) 

where C4 = 4.2 x 10^. The synchrotron self-absorption optical 
depth can be obtained tq = K5/3(Xj,i,)/Ks/3(x^^) « (ly^-Jiypk)^^^, 
where z^sa is the self-absorption frequency. Following Ghis- 
ellini et al. (1993), we take the synchrotron self-absorption 
Vsa = 3.0GHz in BL Lacs. The radiation intensity, in the co- 
moving frame, is given by, = nQCT,BT(x)RQ, which is as- 
sumed to be homogeneous and isotropic. The energy density 
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of radiation field can be written as f/^^j = 4ttc ' I^idv' k, 
ATTC^^c^noBRov'^, where we have used F{x) « S{x-xo) and 
xo = 0.29. The avoidance of Compton catastrophe leads to a 
condition 

uL = xM, (4) 

where Xi < 1' ™d the magnetic field energy density, = 
B^/Stt. In fact the condition Xj = 1 is generally satisfied in 
most BL Lacs. Lastly, we assume parameter X2 represents a 
departure from equipartition of energy between the magnetic 
field and electrons, 

no-fonieC^ = XiUb- (5) 

X2 = 1 means the equipartition between electrons and magnetic 
field. We are not sure the value of Xj, but we can constrain it 
from the general consideration of unified scheme of BL Lacs. 
Combining the above equations we have the Doppler factor as 



r)_ 19 9, 13/48^5/32 ^1/4 



(6) 



where fpk,i5 = fpk/lO^^Hz, Lpk,44 = Lpk/lO^^erg s \ S = 
^-17/8^5/4^5/12^ ^^^^^ ^ r/sa/3GHz, and typically 2?=12 for 
S = 1. The uncertain parameters Xi and Xi are absorbed into S 
as well as Vsa, which may be different in HBLs and LBLs. Ac- 
cording to the unified scheme of BL Lacs (Celotti et al. 1993, 
Urry & Padovani 1995), typically, the mean Doppler factor in 
HBLs is (I>h) = 6.0 and (Pl) = 10 in LBLs. From this argu- 
ment, we obtain (Sh) = 2.2 x 10""* and (Sl) = 0.74 from the 
sample of Costamante et al. (2001). If x, = 1, namely the pres- 
ence of equipartition between magnetic and radiation fields, we 
have (x^^) = 1.2 x 10"-^ and (x^) = 0.8. This strongly implies 
that the energy density of electrons in HBLs is far away from 
equipartition with magnetic field while this equipartition gen- 
erally holds in LBLs. This may reflect the reason why the peak 
frequency in HBLs are more highly variable than that in LBLs 
(Costamante et al. 2001), because all likely mechanisms of 
electron acceleration relate to magnetic field. We find the di- 
mension Ro oc t'pk'^'^^ipk''^^xy^^x7^^'^> and the number density 

of electrons no oc i^II^^lJ^^^^^x'^'^^^^X^^^- 

The ratio of kinetic to radiative (in comoving frame) lu- 
minosity, Er, can be expressed by = Lkin/^^rad - T^'^J-^Hn/L-p^, 
where L[.^^ is the radiative luminosity in jet comoving frame 
and Lkin = irRlnampC^ (iV'^ (Celotti & Fabian 1993), here fi, Y 
and Mp are the velocity, Lorentz factor and proton mass, respec- 
tively (we assume jet is proton-electron plasma). The number 
density of electrons, dimension of region and Doppler factor 
are already expressed by the peak frequency and luminosity, 



we have e, oc v ! L 



1/4,9/32 107/32 25/16 

A, 



pk %k X 




We then get 



f{L\)y^ 25/16 



ixk) 



(7) 



^pk/ / \ \^pk/ 

where we take Xi = 1- We have (e^) = 1.5 x 10"'*(e^), which is 
obtained from the sample of Costamante et al. (2001). Inter- 
estingly the derived ratio of kinetic and radiative luminosity is 
roughly consistent with the results in Celotti & Fabian (1993), 
showing a bimodal distribution of this ratio. 

It was originally suggested that an optically thin torus powers 
the central engine of BL Lac (Rees et al. 1982), as stressed re- 
cently by Cavaliere & D'Elia (2001) in their model of blazar 
main sequence. The powerful emission from the relativistic 
jet originates from the efficient extraction of energy from the 



black hole spin by the Blandford-Zanjek process. Although it 
has been shown that this process in geometrically thin disk is 
not so efficient to pump the energy from black hole spin (Livio, 
Ogilvie & Pringle 1999), it may dominate over the emission 
from disk in the ADAF regime. With the self-similar solution 
of the optically thin ADAF of Narayan & Yi (1994), the output 
power of the BZ process is 

Lbz = eBzMc^ with ebz « 1.8 x 10"^ (8) 
where M is the accretion rate of a black hole with maximum 
spin (Armitage & Natarajan 1999). The ADAF dominates when 
the accretion rate is below the critical value < o^Mem, 
where a is viscosity parameter and the Eddington accretion rate 
is MEdd = A-ttGMbh/ VKenC = 1 .64 X 10^''MBH,8(g/s), the electron 
scattering opacity Kes = 0.34, = 0.1 andMBH.s =Mbh/1O^M0 
(Narayan & Yi 1995). This condition works within a few hun- 
dred Schwartzchild radii for massive black holes, for which it is 
independent of the black hole mass (Narayan & Yi 1995). Ac- 
cording to the model of their main sequence, the blazars' phe- 
nomenon can be understood in terms of different levels of m. 
The accretion rate m in BL Lac should be of 10"^^"""*' (Cavaliere 
& D'Elia 2001) (here we use different definition of Eddington 
limit). For all the BL Lacs, it may be a good approximation that 
the critical value of the dimensionless accretion rate may be the 
same ;w « a^. This condition then yields a lower limit of black 
hole mass via L'^^ = Lbz/ Cr. 

M^r>i ''l^'^i^^S^^''- (HBLS), 

-\3.22xlOU^;^/IX£,eMe (LBLs), 
where we have used Lpk = f^i^a (e.g. Sikora et al. 1997), and 
^ = (eL)5/a^2. Here Q!_2 = a/0.01, and (ef')5 = (ef')/2.5 x 10^. 
a = 0.01 is reasonable for supporting ion tori (Rees et al. 1982) 
since a high a. ADAF does not possess empty funnels (Narayan 
et al. 1998). It should be addressed that the equation (9) is 
mainly based on the blazar main sequence suggested by Cava- 
liere & D'Elia (2001), namely, m^a^= lO-'^^"*', which should 
be satisfied for BL Lac objects as a special group of blazar ac- 
cording to the accretion level. 

On the other hand, Bondi accretion rate sets an upper limit of 
supermassive black hole mass. Chandra's observation shows 
that three giant elliptical galaxies, NGC 1399, NGC 4472 and 
GNC 4636, have far below Bondi accretion rate using the ob- 
served hot interstellar gas surface brightness profiles (Loewen- 
stein et al. 2001). This strongly hints that the activity of 
giant elliptical galaxies may be due to the level of accretion 
rate, namely the activity condition is m > meondi- The dimen- 
sionless accretion rate given by equation (9) is m = 6.88 x 

l^ '^'^pk^ie'^^pk^^iH.V' where we take ^ = 1. The Bondi ac- 
cretion rate in giant elliptical galaxies is 

WBondi = 1.1 X lO^'^MBH.gn-ic;], (10) 

where ri-i = n/O.lcm"^ is the number density of the interstellar 
medium and Cs,5 = Cs/500Km s"^ is the sound speed (Di Matteo 
& Fabian 1997). The upper hnoit of black hole mass is given by 
the condition rh > WBondi. 

Mir < 1-5 X I0\^r Lli%nTc%' M^. (11) 
The black hole mass is then confined by 

M^g^P<MBH<M^™<". (12) 
The physical meanings of equations (9), (11) and (12) reflect 
the intrinsic relation among the central engine, the emission 
and even its host galaxy, which aUows us to deduce the mass 
of black hole. 
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3. MASSES OF BLACK HOLES IN BL LACS 

There is a largest sample including the high and low 
frequency-peaked BL Lacs composed of 1 Jy and Slew Sur- 
vey BL Lacs, in which the peak frequency and its luminosity 
are given by fitting the SED with a polynomial (Costamante et 
al. 2001). There are three exceptional objects, Mrk 501, lES 
2344H-514 and lES 1426H-428 that show highly variable prop- 
erties. Their peak frequencies change by as much as lO'^Hz to 
lO'^Hz and peak luminosities from lO"*^'^ to lO'*"* -' erg/s. We 
exclude the data of these three objects in their high states to 
leave 100 data. All the data can be obtained from Costamante 
(private communication). 

Figure 1 shows the distributions of upper and lower limit of 
black hole masses in BL Lac objects. The uncertainty in eq.(9) 
is due to the parameter ^, namely, a and (e^). However we 
can use their typical values, a = 0.01 and (ejr) = 2.5 x 10^, then 
^ = 1 . We find the high mass peak is in good agreement with 
Laor's limit, which he find the mass of black hole in radio-loud 
quasar should exceed lO^M© (Laor 2000). Thus we think ^ = 1 
is reasonable. 
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Fig. 1 . — The mass distribution of blaclc holes in BL Lac objects (solid line). 
The high (i^p^ > lO'^Hz) and low (Upi^ < lO'^^Hz) frequency-peaked BL Lacs 
ai'e usually refeired to as blue and red BL Lacs, respectively. The upper pannel 
is the lower limit mass distribution whereas the low pannel is the upper limit 
mass distribution. There is a bimodal mass population of black holes in HBLs 
and LBLs. It is very interesting that the mass distribution of black holes in nar- 
row line Seyfert 1 galaxies (dotted line) overlaps that of high frequency-peaked 
BL Lac objects. 

It is interesting to find there is evidently a low peak of black 
hole mass populations centered roug hly at 7.0 X lO^M© from 
(1.0 30) X 1O^M0 when the high peak reconcile with the 
Laor's limit. The mass of black hole in Mrk 421 has been es- 
timated to be w 1O^M0 with accretion rate m ^ 10~-'(Celotti, 
Fabian & Rees 1998). The present result overlaps this mass. 

Many intermediate BL Lacs have been found in the RGB 
(Laurent-Muehleiser et al. 1998), DXRBS (Perlman et al. 
1998) and REX surveys (Caccianiga et al. 1999). Simulations 



show that the IBL content, even the deeper ones, is affected 
more strongly by the sample flux than by intrinsic properties 
(Giommi et al. 2000). We have not an homogeneous complete 
sample composed of HBLs, IBLs and LBLs. It is still unknown 
the unbiased occurring of BL Lacs according to peak frequency 
Upi,. It may have two possible cases. First if the IBL fraction to 
the total of BL Lacs is still smaller than that of HBL and LBL, 
the bimodal distribution of BL Lacs still holds. Second, this 
fraction is comparable with the other two, then the gap between 
the two peaks will be filled up. There is no available data of 
peak frequency and peak luminosity for a sample of interme- 
diate BL Lacs, but we presume an arbitrary distribution of the 
second case of IBL for illustration in Figure 1 . 

4. BLACK HOLES AND SPHEROIDS IN BL LAC OBJECTS 

The black hole mass is generally thought to increase with 
the mass of spheroid in galaxies (Magorrian et al 1998) and 
quasars (Laor 1998), especially the mass is strongly correlated 
with dispersion velocity (Ferrarese & Merritt 2000, Gebhardt et 
al 2000). The high mass family of black holes in LBLs roughly 
agrees with the Magorrian et al. relation. 

With the mean absolute magnitude of X-ray selected BL 
Lacs at R band (Mr) = -23.9 (Falomo et al. 2000), we get 
the mean mass of black holes in HBL (M™^) = IO'^Mq if 
the relation log(MBH/MQ) = -0.51Mr- 3.2 (McLure & Dun- 
lop 2001) works. Clearly, the low mass family consisting of 
HBLs does not satisfy the Magorrian et al. relation. If the mass 
of black hole is of lO'^M© in HBLs, the dimensionless accre- 



tion rate from equation (9) is m = 6.88 x ^0~^v^l\^^^L^^^sM^li^g, 
We should note that: 1) Such a low accretion rate breaks down 
the blazar's main sequence according to accretion level (Cav- 
aliere & D'Elia 2001). 2) We find m < jiiBondi (from eq.lO) 
if the mean mass of black holes in HBLs is (Mbh) = IO^Mq. 
If the accretion rate in HBLs is indeed far below Bondi accre- 
tion rate, this means that all the face-on giant elliptical galax- 
ies will appears as HBLs. If the maximum viewing angle 
is 6'max, the appearance frequence of BL Lac objects will be 
given by / = ( 1 - cos 6'„ax)/2 ~ 6.7% with 9^^^ = 30° (Urry & 
Padovani 1995) based on random orientation of giant ellipti- 
cal galaxies. The number density of BL Lacs will be given by 
A^BL = fNcE, where A^ge is the elliptical galaxy number den- 
sity. We get the number density of giant elliptical galaxies is 
conservatively A^ge ^ 10"* Gpc"^ for My < -23.0 in the local 
space (e.g. Table 3 in Woltjer 1990). The predicted number 
density of BL Lacs based on the giant elliptical galaxies will be 
fNcE ^ 600 Gpc~^. However, the observed number density of 
BL Lacs is only A'bl < 1 Gpc"^ estimated from the luminosity 
function in their Fig. 18 of Urry & Padovani (1995). Evidently 
the predicted number density is much higher than the observed 
by a factor of > 600! To iron out this inconsistency is to lower 
the mass of black holes in HBLs. These arguments suggest that 
the black hole mass might not obey the Magorrian et al. re- 
lation and support the blazar main sequence according to the 
accretion level. 

There several narrow line Seyfert 1 galaxies as radio-loud ob- 
jects, which ai-e PKS 2004-447 {R > 1700) of 5x lO^'M© (Osh- 
lack, Webster & Whiting 2001), PKS 0558-504 (R = 27), RGB 
J0044H-193 {R = 31) (Siebert et al. 1999) and RX J0134.2-4258 
(R = ll) (Grupe et al. 2000). These objects evidently comprise 
a special population in the plot of radio loudness and black hole 
mass (Laor 2000). The distinguished features of narrow line 
Seyfert 1 galaxies (NLSIs) can be explained by the hypothesis 



4 



Wangetal. 2001 



that they contain less massive black holes with relatively high 
accretion rates, close to or even super-Eddington, and disks 
which are face-on to the observer (Boiler, Brandt & Fink 1996) 
as BL Lacs. A heterogeneous sample of 59 NLSl galaxies 
with good spectroscope observations is provided (Veron-Cetty, 
Veron & Goncalves 200 1 ) . The masses of black holes in NLS 1 s 
has be obtained from their H/3 width (Wang & Lu 2001), we 
plot the mass distribution in Figure 1 . It is found that the mass 
distribution of black holes in high frequency-peaked BL Lacs 
nicely overlaps that of the narrow line Seyfert 1 galaxies. We 
thus suggest that there might be an evolutionary relation from 
a BL Lac to a narrow Une Seyfert 1 nucleus although their host 
galaxies are quite different. 

5. CONCLUSIONS AND DISCUSSIONS 

This study shows that there are two populations of black 
holes in BL Lacs corresponding to the low mass family of HBLs 
and the high mass family of LBLs. The high mass family sat- 
isfies the Magorrian et al. relation roughly and is in agreement 
with Laor's limit whereas the low mass family may not obey 
the Magorrian et al relation. 

Based on the unified scheme of HBLs and LBLs, we find 
that the energy density of electrons is far below the equipar- 
tition with magnetic field in HBLs whereas this equipartition 
roughly holds in LBLs. This reflects why the peak frequencies 
in HBLs are usually highly variable. 

The bimodal distribution of black holes sets strong con- 
straints on the model of evolution of BL Lacs (Ghisellini et al. 
1998, Bottcher & Dermer 2001). First if accretion rate of black 
holes, rii » Q^, holds in HBLs constantly, it will take more than 
lO^yrs for black holes to grow up to be the same as in LBLs. 
In such a case it seems that there is no evolutionary connection 
between HBLs and LBLs. On the other hand, if the model of 
multiple generation of short-Uved quasars works (Haehnelt et 
al 1998, Kuhn et al 2001), the accretion rate increases and de- 



creases around its peak (Cavaliere & Vittorini 1998), and then 
repeats in the life. HBLs are only short-lived objects, there 
might be evolutionary sequence from HBLs— >LBLs. The pos- 
sible intermediate objects are narrow line Seyfert 1 galaxies, 
which are undergoing very fast accretion. The mass distribu- 
tion of black holes in high frequency-peaked BL Lac objects is 
similar to that of the narrow line Seyfert 1 galaxies, likely im- 
plying some evolutionary cormections between the two kinds of 
galactic nuclei. 

The most uncertainty is the parameter ^ due to a and e^, but 
the present results can not be modified since ^ can not changed 
three orders. On the other hand ^ can be constrained by the 
Laor's limit, which makes the low mass family safe. It is gen- 
erally thought the relativistic jet is originated in the vicinity of 
central black hole. We would like to stress here that the emis- 
sion from jet may directly take the most important information 
of black hole parameters. The present method is robust in the 
view of this point. Although we use the monoenergetic distri- 
bution of relativistic electrons to calculate the emission from 
jet, the present method works as a robust tool to extract statistic 
properties of BL Lacs. It highly desired to measure the mass 
of black holes via stellar dynamics or gas dynamical methods 
in centers of BL Lacs to confirm the present results. The orig- 
ination of the bimodal distribution of black hole masses in BL 
Lacs remains open. 
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